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ABSTRACT 

We derive stellar masses from SED fitting to rest-frame optical and UV fluxes for 401 star-forming 
galaxies at z ~ 4, 5, and 6 from Hubble- WFC3/IR observations of the ERS combined with the 
deep GOODS-S Spitzer/IRAC data (and include a previously-published z ~ 7 sample). A mass- 
luminosity relation with strongly luminosity-dependent M /~Luv ratios is found for the largest sample 

(299 galaxies) at z ~ 4. The relation A4 oc Lyy has a well-determined intrinsic sample variance 
of 0.5 dex. This relation is also consistent with the more limited samples at z ~ 5 — 7. This z ~ 4 
mass-luminosity relation, and the well-established faint UV luminosity functions at z ~ 4 — 7, are 
used to derive galaxy mass functions (MF) to masses Ai ~ 10 8 at z ~ 4 — 7. A bootstap approach 
is used to derive the MFs to account for the large scatter in the Ai-Luv relation and the luminosity 
function uncertainties, along with an analytical crosscheck. The MFs are also corrected for the effects 
of incompleteness. The incompleteness-corrected MFs are steeper than previously found, with slopes 
olu ~ —1.4 to —1.6 at low masses. These slopes are, however, still substantially flatter than the MFs 
obtained from recent hydrodynamical simulations. We use these MFs to estimate the stellar mass 
density (SMD) of the universe to a fixed Mtjv.ab < —18 as a function of redshift and find a SMD 
growth cx (l+z) _3 - 4±0 - 8 from z ~ 7 to z ~ 4. We also derive the SMD from the completeness-corrected 
MFs to a mass limit M. ~ 10 8 M Q . Such completeness-corrected MFs and the derived SMDs will be 
particularly important for model comparisons as future MFs reach to lower masses. 

Subject headings: galaxies: evolution — galaxies: high-redshift 

1. INTRODUCTION 

Measurements of the stellar mass of high-redshift 
galaxies (z > 4) provide important constraints on sce- 
narios of galaxy formation and early evolution. For 
example, the evidence for a strong correlation between 
the observed star f ormation rate (SFR ) and stellar mass 
(jStarket al.ll2009t lLabbe et alJl2010al lbl). wi th little ap- 
parent evolution between z ~ 3 and z ~ 7 (jStark et al.l 
2009; Gonzale z et al.|[20i0f ). suggests an epoch of expo- 
nential growth which, in terestingly, is simila r to that 
found in simulations (e.g. iFinlator et aTll2010[ ). 

Recent d eep near-IR WFC3/IR observations over the 
ERS field (jWindhorst et al.l l2010t ) combined with pre- 
existing deep GOODS IRAC data provide access to the 
rest-frame UV and optical wavelengths of 4 < z < 7 
star forming galaxies and hence reasonably accurate es- 
timates of their A4/L ratios and stellar masses. The 
substantial samples of z > 4 galaxies detected with 
WFC3/IR span a range in stellar mass, allowing the 
derivation of mass functions (MFs). 

MFs are fundamental characteristics of the galaxy pop- 
ulation but in practice, they are difficult to compute di- 
rectly, especially at high redshift because of selection ef- 
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fects, incompleteness, and contamination by interlopers. 
An alternative approach to derive the MF is to start with 
the well-determined UV LF at these redshi fts and convert 
to ste llar mass using an average Ai/L (e.g. lMcLure et all 
2009). The main advantage is that LFs are corrected for 
all selection effects in the data and reach very faint lim- 
its, although such a simple approach does not take into 
account any possible luminosity dependence or scatter 
of the A4/L. In this Letter, we estimate improved MFs 
at 4 < z < 7 from the published UV-LFs by deriving a 
relation between UV luminosity and stellar masses. We 
also estimate the scatter in the Ai/L, which allows us to 
correct for incompleteness at low stellar masses. We use 
the MFs to compute the stellar mass density (SMD) of 
the universe at z = 4, 5, 6 and 7. 

We adopt an CIm = 0.3, f^A = 0.7 cosmology with 
H = 70 kms s^ 1 Mp c" 1 throughout. A ll magnitudes 
are in the AB system (|Oke fc GunnlH98l . 

2. GALAXY SAMPLE FROM HST AND SPITZER DATA 

The sources used here for determinations of the Ai/L 
ratios and z — 4 — 6 MFs were found in the recent 
Hubble- WFC3/IR observations of the ERS field. Both 
the GOODS ACS optical (-B435 ^06*775^850) an d the 
WFC3/IR (V g 8 Jno-ffi6o) data rea ch depths of ~ 28 mag 
(5a, 0.35"-diameter a pertures, see lBouwens et all 120101 : 
IGiavalisco et all 12004). All sources have Spitzer/IRAC 
coverage with depths of 27.8 and 27.1 in the [3.6] and 
[4.5] channels, respectively (la in 2.4"aperture s). The 
z ~ 7 sample is taken from ( Labbe et alll2010aj ). 



The z ~ 4, 5, 6 sample totals 679 objects, consisting 
of 524 B, 123 V, and 32 i-dropouts that were selected with 
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Fig. 1. — Stellar masses as a function of UV Luminosity (Mtjv.isoo = 51.63 — 2.5 X log 1 Q(Luv,i50o[ er g s Hz ])) for the z ~ 4 , 5, and 
6 samples. SFR uncorr (top axis) is derived using the lMadau et al.l JT998I') conversion formula (no extinction correction). The final sample of 
401 sources with FAST SED-fit mass estimates is shown here. Open squares indicate low-S/N measurements (< 2a in [3.6]). The darker 
symbols in each panel represent the median mass of the sample (~ 0.5 Muv,i500 _ma g bins). The small error bars represent the bootstrapped 
errors. The larger black error bars include a conservative estimate of the systematics computed by comparing the estimated median mass 
at a given luminosity with the mass estimated from the stacked SEDs at the same luminosity. The dashed blue line (slope = 1.7 ± 0.2) 
represents the mean trend between logioA4 and Mtjv 1500 at z ~ 4. It is consistent with no evolution with redshift. The scatter at the 
luminous end (±0.5 dex), where photometric errors are small, is intrinsic (see Figure l2ll. 



the same criteria as iBouwens et al.l (|2007ft : 
z ~ 4 B-dropouts: 

(B435 - ^606 > 1-1) A [B435 ~ ^ 6 06 > (V606 - 2 85 o) + 1-1] 
A {V 606 " Z850 < 1.6) 

z ~ 5 V-dropouts: 

{[^606 - «775 > 0.9(i 775 - Z 850 )) V (V 60 6 - 1775 > 2)} 

A (V 6 oa - i 775 > 1.2) A (i 775 - z 850 < 1.3) 
z ~ 6 i-dropoutf0: 

(«775 - 2850 > 1-3) A (z 8 50 ~ ^125 < 0.8) 

The rest- frame optical photometry from Spitzer/IRAC 
is ideally suite d for deriving stellar masses at these 



redshifts (e . g. iPapovich et alj [20011: lYan et al.l [20051: 
lLabbe et al.l 12010a ) . A challenge is that the broad 



IRAC PSF usually results in these faint sources be- 
ing contaminated by foreground neighbors. To obtain 

6 Slightly modified, Bouwens et al. 2010, in prep. 



rel iable IRAC fluxes we use th e deblending method 
of lLabbe et all (2 006. see also i Gonzalez et al.1 120101: 
lLabbe et al. 2010b]i; IWuvts et al. |[2QQ7|;|de Santis et al.l 
|2007|) . Briefly, this method uses the higher- resolution 
HST images to create models of both the foreground 
neighbors and the source itself. We convolve each model 
image with a kernel to simulate the IRAC observations. 
We fit for all the sources simultaneously (with indepen- 
dent normalization factors) and subtract the best fits for 
the neighbors. In the clean image of each dropout we are 
able to perform standard aperture photometry. We use 
2.5"-diameter apertures and correct the fluxes to total 
assuming stellar profiles (1.8 x in both channels). 

As expected, our cleaning procedure does not work 
for every source. We restrict our sample to the 60% of 
sources with the best x 2 residuals. This reduces the num- 
ber of non-optimal subtractions to < 8%. The final sam- 
ple suitable for deriving masses from the HST+Spitzer 
data totals 401 sources: 299 z ~ 4, 78 z ~ 5, and 24 at 
z ~ 6. We do not expect this selection step to introduce 
any important biases, since it depends on the distribu- 
tion of the non-associated neighbors of the source. Of the 
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remaining sources, ~ 50% have low IRAC S/N (< 2 a in 
[3.6]). 

3. STELLAR MASS ESTIMATES FROM SED FITS 

We use the FAST SED-fitting code (|Kriek et al.ll200ll 
to derive stellar masses of the 401 z ~ 4 — 6 sources from 
SED fits with the full suite of fitted parameters. For 
all sources we fit the broadban d ACS+WFC3/IR+IRAC 
[3.6] and [4.5] fluxes u sing thelBruzual fc Chariot! (pOOl 
BC03) models with a lSalpeterl (1 19551 " initial mass func- 
tion (IMF, 0.1-100 M ) and assuming a 0.2 Z Q metallic- 
ity. We also include the sample of z ~ 7 galaxies wi th 
similarly-determined masses from lLabbe et all (|2010af >. 

The SFH cannot be uniquely determined from broad- 
band SEDs due to well-known degeneracies between the 
star formation timescale, age, and dust extinction. We 
have assumed a SFH with a constant SFR. Different 
SFHs introduce systematic offsets to the mass determina- 
tions , largely independent of redshift (cf. iPapovich et al.l 
I2010I ). The systematic differences between masses based 
on decli ning, constant, or r ising SFHs are typically < 
0.3 dex (jFinlator et al.ll2007l ). 

Figure [T] (left) shows the FAST SED-fit stellar masses 
(from HST+Spitzer data) versus UV luminosity (bottom 
axis). While the scatter is large (RMS^ 0.5 dex), there is 
a clear trend of increasing mass with increasing UV lumi- 
nosity. The stellar mass M.-\^jjv,i^m relation at z ~ 4 is 
well-fit by logio(X) oc 1.7(±0.2)'logio(Li 5 oo). The lower 
bound that appears at masses < 10 s M Q corresponds 
to the Ai/L of the youngest model we allow (lOMyrs). 
This is a reasonable assumption for the majority of the 
sample, as IRAC detections in the stacks suggest that 
continuum is dominant, not emission lines. This is also 
not critical, since our inferred Ai/L trend is insensitive 
to the cutoff. The z ~ 4 relation is consistent with the 
z ~ 5 sample, and, in zero-point, with the small z ~ 6 
sample, and also w ith the z ~ 7 sample presented in 
(jLabbe et al.ll2010aft . 

Figure [2] explores the z ~ 4 Ai/L ratio trend in more 
detail, showing that Ai /Lyv',1500 depends on luminosity; 
the Ai/L ratio is ~ 5x lower at M{/y,i500 = —18 than at 
M[/v,i500 = — 21. This suggests that UV-faint galaxies 
contribute l ess to the global SMD than assumed in previ- 
ous studies ()Labbe et al.ll2010"bU af). Howe ver, due to the 
steep faint-end slope of the UV LF (e.g. iBouwens et al.l 
I2010T ) . their high number density makes their total con- 
tribution quite significant. 

A striking aspect of the relation is the large scat- 
ter in Ai/L. The observed sample variance (one stan- 
dard deviation) for our sample is ~ 0.5 dex for 
-21 < Muv,i500 < -18- At the bright end M UVt i 500 < 
—20 the scatter is largely intrinsic, whereas at the faint 
end Myy^isoo > —19.5 it is dominated by observational 
uncertainties. In particular, the stellar masses of sources 
with IRAC detections are much better constrained than 
IRAC-undetected sources. Photometric uncertainties 
contribute ~ 0.14 dex to the scatter at M[/v,15qo ~ —20 
(0.37 dex at —19). Moreover, we find that the Ai/L ra- 
tio is tightly correlated with the J — [3.6] color (standard 
deviation 0.18 dex, Figure [21 right), suggesting that the 
variation is real, and not an artifact of the modeling. 

The relation in Figure [5] (right) also allows us to es- 
timate the possible effect of contamination by emission 
lines (not included in our models). At z ~ 4, a 20% con- 



tribution of Ha to [3.6] would result in redder J — [3.6] 
colors and hence overestimates of the Ai/L and of the 
masses by 30%. This would affect the SMDs at all red- 
shifts because they all rely on our z ~ 4 Ai/L ratio 
estimates (see §4). 

4. STELLAR MASS FUNCTIONS ATZ~4, 5, 6, AND 7. 

The stellar MFs at a given redshift can be estimated 
by combining the LF at that redshift with an appropri- 
ate Ai /L relation. Since UV LFs have been derived from 
large samples to very faint limits and carefully corrected 
for a wide range of potential biases, they constitute an 
excellent basis for determining MFs. The calibration 
of log(A'l) vs M[/y,i5oo in Figures Q] and [2] provides the 
means to transform the z ~ 4 — 7 UV-LFs into MFs. 

However, the scatter about the mean A't-Mj/^isoo re- 
lation is so large that ignoring it would produce signifi- 
cant errors. Galaxies with relatively low luminosity but 
high Ai/L ratios, for example, contribute significantly at 
the high mass end of the MFs. Hence, we take care to 
determine the average A4-M;7y,i500 relation in a robust 
way, we characterize the scatter at the high mass end 
and use this estimate of the scatter at lower luminosi- 
ties/masses where the observational uncertainties domi- 
nate. 

We use two approaches to create the MFs. First, we 
use the individual points in Figure [TJ as representative 
of the A'f-Mf/^isoo distribution by bootstrap resampling 
them. To correct for incompleteness we add faint sources 
to the distribution. This important step increases the low 
mass slope of the MFs substantially. Second, we use the 
fitted A'l-Mj/v',1500 relation and an idealized model of its 
scatter to produce what we label as "analytic" versions 
of the MFs. We compare to other estimates as a cross- 
check. 

Boo tstrapped MFs: We start wi th the z ~ 4- 7, UV- 
LFs of IBouwens etail (|2007l [2010h and draw 40000 lumi- 
nosities from each LF in the range —21.5 < Muv,i500 < 
— 18. We convert the luminosities to stellar masses by 
bootstrap re-sampling from the distribution of points at 
z ~ 4 in Figured] We use the z ~ 4 log(Af) - M (7^,1500 
distribution at all redshifts because it is well-defined over 
a wide range of luminosities and is consistent with the 
relations at other redshifts, including the z ~ 7 rela- 
tion presented in lLabbe et al.l (|2010a[) . To account for 
the uncertainties in the LFs we perturb their Schechter 
parametrizations within the uncertainties and repeat 
5000 times. This "bootstrap" process results in the un- 
corrected MFs (Figure |3] open squares and error bars) . 

As a crosscheck we also derived a histogram MF at 
z ~ 4 directly from the masses of the z ~ 4 sample using 
the search volume for the i?-dropouts. This straightfor- 
ward process gives a MF that is identified in Figure [3] 
as the "direct MF" for comparison with the uncorrected 
"bootstrap MF" . 

To correct the MFs for incompleteness at Ai < 
1O 8 5 M0, we re-derive the MFs but now including fainter 
(—18 < Muv,i5oo < —15) sources. We extrapolate the 
observed log(A / f)-M[/y : i 5 oo relation to lower luminosities 
and assume that the low-luminosity scatter is similar to 
the scatter around Muv,i500 ~ —18.5. The resulting er- 
rors on the corrected points include an added uncertainty 
that is typically about 30-40%. This accounts for the 
LFs uncertainties (specially the faint end slope) and the 
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Fig. 2. — (left) -M/L ratio as a function of UV Luminosity for the z ~ 4 sample. Symbols and error bars as in Figurc[T] The median A4/L 
ratio changes by a factor 5x in the luminosity range our sample, (right) The correlation between the M/L and J — [3.6] color. Arrows 
indicate 2 a- upper limits. This tight relation suggests that the large scatter observed in the A4/L (left panel) is largely due to intrinsic 
variations in the UV-to-optical colors. Photometric scatter can only account for < 0.14 dex at Mtjv 1500 ~ — 20 (0.37 dex at —19). 



large scatter about the relations. Other sources of un- 
certainty may remain but further assessment is needed 
to fully evaluate them. Regardless, the current correc- 
tions must make these corrected MFs a better estimate 
of the true MFs. Applying the completeness corrections 
is a crucial step and significantly changes the slope of 
the MF at lower mass. The corrected MFs are shown in 
Figure [3] by the solid points and the solid color band, and 
are referred to as the bootstrap MFs. 
Analytic MFs: The best fit M.-Muv,i500 relation in 

Figure [I] at z ~ 4 is: M oc L^y 1500 ■ We combine 
this relation with the same LFs as above to derive the 
analytic MFs (Figure [3]) . The large scatter needs to be 
accounted for in order to generate a realistic estimate. 
The mass distribution at a given luminosity is assumed 
to be log-normal with a standard deviation of 0.5 dex, 
centered on the fitted \og(M)-M.uv,isoa relation. The 
resulting relation, normalized by the LF, is integrated 
over L to get the analytic MFs. 

At lower masses, the slope of the MF is set by the faint- 
end slope of the UV-LFs and of the log(A / /)-Mf/y i i 50 o 
relation. Extrapolating this M-M-uv.isqq relation to 
lower luminosities results in steep low-mass slopes for 
the MFs of -1.43(±0.11), -1.39(±0.11), -1.44(±0.15), 
and -1.55(±0.21) at z ~ 4, 5, 6, and 7. These 
MF slopes are sli ghtly flatter than the UV LF slopes 

(a = -1.7 2: iBouwens et "ail 120101 ) . They are in 

good agreement with our completeness-corrected boot- 
strap MFs, and so provide a useful "sanity check" on 
those results. The standard deviation of 0.5 dex in the 



log(A / l)-Mc/y i i5oo relation results in a slightly enhanced 
number density at the high-mass end compared to a case 
with no scatter. Our corrected MFs are considerably 
steeper than othe r MF determinations at high redshift 
(Star k~et al.l 120091 ) that do not apply completeness cor- 
rections. Truncating the analytic M.— Mj/y^oo relation 
at M1500 < — 18 (to represent survey incompleteness) re- 
sults in analytic MFs that are in good agreement with 
the non-corrected bootstrap MFs at low masses. 

Comparisons to simulated MFs: The MFs derived here 
are substantially steeper at low masses than what has 
been found in the past at these redshifts. This is not 
unexpected given that it has not been typical to correct 
for incompleteness. While the corrections are uncertain 
in magnitude the sign of the correction is not. Inter- 
estingly, even with the corrected and steeper slopes, the 
observed MFs are quit e different from what is s een in re- 
cent simulations (e.g. iChoi fc Nagaminel 120101 : Finlator 
et al. in preparation). The simulated MFs are steeper, 
with dramatically more low-mass sources than we find. 



5. STELLAR MASS DENSITY ATZ~4, 5, 



AND 7 



The MFs can be integrated to determine the SMD of 
the universe at high redshift. First, we integrate the 
(uncorrected) bootstrap MFs to determine the SMD at 
z = 4, 5, 6, and 7 to faint-luminosity limits (M1500 < 
-18 - Table [TJ Figure H left). Fitting the SEDs us- 
ing the observed fluxes rather than uppe r limits al- 
lows u s to reach lower limits than those of iStark et al.l 
(2009) at z — 4 — 6. To compare to those results, 
we correct their original M1500 = —20 limit to our 
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Fig. 3. — Stellar mass functions at z ~ 4, 5, 6, and 7 derived from the log(A4) — Mtjv,i500 distribution for the z ~ 4 B-dropouts 
(Figure[TJ, and the Bouwens ct al. (2007, 2010) UV-LFs at z ~ 4 — 7. The points are derived from the "bootstrap" approach (see text). 
Frrors reflect uncertainties in the LF and the ~ 0.5 dex Iff scatter of the JV[— Mpy 1500 relation (Figure[TJ. Completeness-corrected values 
are estimated assuming that the M—bAjjV 1500 relation extends to fainter limits with similar scatter about the extrapolated mean trend 
(M[/^i5oo < —18 uncorrected: open; corrected: filled; dark band is at 1 a around the corrected values). The direct MF at z ~ 4 (thick 
histogram) is in good agreement with the uncorrected MF (see text). For masses > 10 9 ' 5 Mq, the uncorrected z < 7 MFs are in rough 
agreement with the determinations of Stark ct al. (2009) and of IMcLure et al.l 120091 ') at z ~ 6 and M > 10 10 M Q . The thick dashed 
curve in each panel represents the analytic MFs derived from an idealized A4 — Mjjv 1500 relation (see text §4). These MFs have low-mass 

slopes otM ~ —1.4 1.6, slightly flatter than the UV LFs (ct = —1.7 2.0: Bouwens ct al. 201(j). In turn, the assumed symmetric 

scatter of 0.5 dex flattens their slopes at the high-mass end. The z ~ 4 analytical MF is repeate d in t he other panels for comparison (thin 
dashed curve). The dotted and thin solid lines show the simulated MFs from Choi & Nagamine (2010) and Finlator et al. (in preparation). 
Our new results are corrected for incompleteness, yet the difference between our results and the simulations is already substantial by 
hA = 10 9 Mq. The source of the disagreement is unclear. 
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Fig. 4. — (left) Stellar Mass Density as a function of redshift for sources brighter than Mtjv,1500,AB = —18. These SMD values are 
derived by integrating the uncorrected bootstrap MFs in Figure [3] to the faint luminosity limit Mjyyisoo = —18 at z = 4, 5, 6, and 7. 
For comparison, we show the SMD determinations from Stark ct al. ( 2009j) corrected from their original M(jv".i500 = —20 limit to our 
M(/yi5oo = —18 limit (see text). The low-redshift open circles were derived by integrating the Marchcsini ct al. (2009) MFs between 
8.3 < logio(.A/f/M0) < 13 and multiplying by 1.6 to match the Salpctcr IMF. A constant SFH and O.2Z metallicity was assumed to 
derive the masses at z > 4. The effect of a possible 20% correction due to contamination by Ha is shown, as is the effect of using a different 
IMF. The integrated mass growth we derive with cosmic time is well fit by logio(SMD) oc (1 +z)" 3 - 4±0 ' 8 . (right) As for the left panel 
but now to a fixed mass limit > 10 8 M Q . The SMD to a fixed mass limit is compared to the flux limited values from the left panel. The 
differences are relatively small (see text). However, the importance of being able to utilize the completeness-corrected MFs will increase as 
improved, deeper data becomes available and we can push to lower mass limits. 
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Summary of Results. 
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M1500 = -18 limit by adding 0.18, 0.22, and 0.32 dex 
at z ~ 4, 5, and 6 respectively. We derive an integrated 
mass growth across cosmic time that is well fit by the 
function logio(SMD) oc (1 + z)- 3 - 4±0 - 8 . The effect of a 
different IMF and of a potential contamination by 20% 
H a at z ~ 5 — 6 is also shown on Figure |U 

A major result of this paper is the derivation of MFs 
corrected for incompleteness at low masses. To utilize 
the (more representative and accurate) corrected MFs 
in deriving the SMD of the universe at high redshift we 
need to integrate to a fixed mass limit. We choose M. > 
10 8 M Q to extend to the limit of our corrected data. The 
right panel of Figure U] shows these results, and compares 
them with the estimates to a fixed luminosity limit (see 
also Table [J). 

The differences are relatively small, less than we ex- 
pected given the size of the corrections. A number 
of effects impact the comparison. For the flux-limited 
(Mab < —18) sample, we include relatively bright 
sources that have masses > 10 s M©. In the mass- limited 
case, we remove those galaxies, and instead include faint 
sources with large masses. Depending on the limits 
adopted and the shape of the LF, this effect is more or 
less important. At z ~ 7, the fainter L* and the steep 
faint end slope make the difference between what is re- 
moved and what is added larger. For the brighter L* at 
lower redshifts, however, the difference is almost null. 

Nonetheless, the derivation of mass-limited SMDs that 
appropriately include corrections at low masses is pre- 
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ferred and will become increasingly important as we push 
our measurements of the MF to lower limits. 

6. KEY RESULTS 

We derive stellar masses from SED fits to HST+Spitzer 
data for over 400 z ~ 4 — 7 galaxies. We determine 
the M-Luv relation and find it to be steep (log(.M) oc 
1.7(±0.2)log(L f/y )) with large intrinsic scatter; the sam- 
ple variance is ~ 0.5 dex at the bright end. We derive 
mass functions by combining the M /L results with pub- 
lished deep UV luminosity functions at z ~ 4 — 7, and 
correct them for incompleteness. The corrected mass 
functions are steeper (a ~ —1.4 to —1.6) than found 
previously, but still far less steep than those from recent 
hydrodynamical simulations. The integrated stellar mass 
density of the universe is derived at z ~ 4, 5, 6, and 7 
to M ~ 10 8 M o . 
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